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A method for the calculation of the electronic structure of pyramidal self-assembled 
InAs/GaAs quantum dots is presented. The method is based on exploiting the C4 
symmetry of the 8-band k ■ p Hamiltonian with the strain taken into account via the 
continuum mechanical model. The operators representing symmetry group elements 
were represented in the plane wave basis and the group projectors were used to 
find the symmetry adapted basis in which the corresponding Hamiltonian matrix is 
block diagonal with four blocks of approximately equal size. The quantum number of 
total quasi-angular momentum is introduced and the states are classified according 
to its value. Selection rules for interaction with electromagnetic field in the dipole 
approximation are derived. The method was applied to calculate electron and hole 
quasibound states in a periodic array of vertically stacked pyramidal self-assembled 
InAs/GaAs quantum dots for different values of the distance between the dots and 
external axial magnetic field. As the distance between the dots in an array is varied, 
an interesting effect of simultaneous change of ground hole state symmetry, type and 
the sign of miniband effective mass is predicted. This effect is explained in terms 
of the change of biaxial strain. It is also found that the magnetic field splitting of 
Kramer's double degenerate states is most prominent for the first and second excited 
state in the conduction band and that the magnetic field can both separate otherwise 
overlapping minibands and concatenate otherwise nonoverlapping minibands. 



Introduction 

Semiconductor quantum dots made by Stranski-Krastanow growth have attracted great 
interest over the past few years from the view of fundamental physics, as well as due to 
their application in optoelectronic and microelectronic devices. In order to understand the 
physics of quantum dots and model and design such devices the electronic structure needs to 
be accurately known. The large lattice mismatch between InAs and GaAs has enabled the 
fabrication of quantum dots putting it at the forefront of both theoretical and experimental 
research. Different quantum dot shapes (such as pyramid^, lens^ and disk-) of InAs/GaAs 
self-assembled quantum dots are often reported. 

A range of theoretical approaches has been used so far to calculate the energy levels 
in self-assembled quantum dots - effective mass^^, k ■ ^ M il & il and the pseudopotential 
methocr^ 1 ^. In quantum dots with cylindrical symmetry, symmetry considerations have 
been applied to effectively reduce the geometry of the problem from three-dimensional 
to two-dimensional, both in the effective mass and the k ■ p method (within the axial 
approximation)^. The possible symmetries of the states in hexagonal Ill-nitride quan- 
tum dots have recently been determined^. The symmetry of the pyramid has been used in 
the effective mass calculation™, to reduce the size of the corresponding Hamiltonian matrix, 
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however in none of the k ■ p calculations of pyramidal quantum dots has the explicit use of 
the symmetry of the Hamiltonian been reported. The aim of this paper is to exploit the 
symmetry in k -p calculation of the electronic structure of pyramidal InAs/GaAs quantum 
dots. 

The symmetry of a pyramidal InAs/GaAs quantum dot when a full atomistic structure 
is considered is C^v and is lower than the symmetry of the dot's geometrical shape 13 . Due 
to the computational complexity of the pseudopotential methods that take into account the 
atomistic nature of the structure, one often employs the k -p method which is considered to 
be a reliable tool for modeling the electronic structure of quantum dots despite its known 
limitations 8 . The symmetry of the k-p model itself is the symmetry group of the zincblende 
crystal lattice. When the model is applied to pyramidal quantum dots the symmetry group 
is the intersection of the geometrical symmetry of the pyramid shape and the zincblende 
bulk symmetry. Since the pyramid shape symmetry group is a subgroup of the symmetry 
of the zincblende crystal lattice, it follows that the symmetry group of the model is the 
double C4 group^. Two different approaches are used to calculate the strain distribution in 
quantum dots - the continuum mechanical 8 - 1 ^ and the valence force field model 85 - 10 . When 
the strain distribution is incorporated in the k ■ p method, the continuum mechanical model 
preserves the C4 symmetry, while the valence force field model, due to its atomistic nature, 
breaks itr^ -. Nevertheless, the comparisons of the two models have shown that they give 
similar results 8 - 1 ^. In this paper, we apply the 8-band k ■ p method with the strain taken 
into account via the continuum mechanical model. We shall refer to this in the rest of the 
text for brevity as the model. Therefore, the symmetry group of this model is C4. All the 
results presented in Sec. II I II are strictly valid in the framework of such an idealized model. 

In previous years, arrays of vertically stacked quantum dots with 10 or more layers have 
been reportecr^ -. The quantum wire behavior of such structures was theoretically in- 
vestigated in Ref. 0. Most of the theoretical investigations of the InAs quantum dots in a 
magnetic field have focused on the dots with a parabolic confinement potential 2 ^^ 2 ^ 1 ^ 2 ^^ 
or just on the conduction band states^^. More recently, the influence of a magnetic field on 
pyramidal single quantum dots has been investigated in the framework of the k-p method 29 . 

This paper is organized as follows. In Sec. [I] the relations necessary to apply the plane 
wave expansion method 30 to 8-band k ■ p calculation of quantum dots based on materials 
with zincblende crystal symmetry in the presence of strain and external axial magnetic field 
are given. In Sec. [Ill the representations of the symmetry group elements, the reduction of 
the representation to its irreducible constituents and finally the symmetry adapted basis are 
found. The states are classified according to their symmetry and selection rules for absorp- 
tion of electromagnetic radiation in the dipole approximation are found. It is also shown 
how symmetry considerations can be incorporated into the plane wave method. Finally, in 
Sec. [inithe method presented is applied to calculate the electron and hole quasibound states 
in a periodic array of vertically stacked pyramidal self-assembled InAs/GaAs quantum dots 
for different values of the period of the structure and in the presence of an external axial 
magnetic field. 

I. THE PLANE WAVE METHOD 

In the presence of an axial magnetic field the total 8-band k ■ p Hamiltonian is a sum of 
the kinetic part of the Hamiltonian the strain part H B , the modification of the kinetic 
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part due to magnetic field Hb and the Zeeman part Hz 

H = Hk + H s + Hb + Hz- 



(1) 



The state of the system within the framework of the 8-band k-p method is given as a sum of 
slowly varying envelope functions ipi(r) multiplied by the bulk Bloch functions \i) (Eq. H]), 
i.e. 



|*> = 5>(r)|<>. (2) 

i=l 

The eigenvalue problem of the 8-band k ■ p Hamiltonian H can therefore be written as 

J2H ij ^ j (r)=E^ i (r). (3) 



The basis of Bloch functions | J, J z ) (where J is the total angular momentum of the Bloch 
function) used to represent the Hamiltonian is given by~ 
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The kinetic and the strain part of the Hamiltonian in the basis (j3J) are given in Ref. [31 
and the matrix elements of the 8-band k ■ p Hamiltonian in the case of a magnetic field of 
arbitrary direction are given in Ref. |32|. In the case of the axial magnetic field that we are 
interested in, the influence of the magnetic field B can be taken into account by replacing 
ki with ki + f Ai in the kinetic part of the Hamiltonian, where Ai is the i— th component of 

the magnetic vector potential, and by adding the Zeeman term Hz (Ref. l32l ). We take the 
symmetric gauge for the vector potential A = \B (—y, x, 0). 

The plane wave method is based on embedding the quantum dot in a box of sides L x , 
L y and L z (Fig. [1]) and assuming the envelope functions are a linear combination of plane 
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waves 

V^O) = 5]Afc ex P ( ik ' r ) > ( 5 ) 
fc 

with the coefficients A^fc to be determined. The wave vectors taken in a summation are given 
by k = 2tt (g, f^) (m x G {-n x , . . . , nj, 777,3, G {-n y , . . . , n y }, m 2 G {-n 2 , . . . , n z }). 
The number of plane waves taken is thus N = 8(2n x + l)(2n y + l)(2n z + 1), where 2n t + 1 
is the number of plane waves per dimension t (t G z}). Due to the symmetry of the 

pyramid, the embedding box sides L x and L y are taken to be equal (L x = L y ), as well as 
the number of plane waves per dimensions x and y (n x = n y ). 
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FIG. 1: Quantum dot geometry. The width of the base is b, the height h, the wetting layer width 
dwL- The embedding box sides are L x , L y and L z . The center of the pyramid base is taken as the 
origin of the coordinate system. 

After multiplying equation from the left by Iv d 3 r exp (— iq ■ r) , where the inte- 
gration goes over the volume of the embedding box, using (jSJ) and the identity 



/ d 3 r exp 

Jv 



2m 
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one arrives at 



where 



hij(q, k)A j:k = EA itQ , 

j,k 



hij(q, k) — — / d 3 r exp (— iq ■ r) Hij exp (ik ■ 



(6) 
(7) 
(8) 



The material parameters in Hk, H s , Hb and Hz are all spatially varying functions since 
they have different values in the quantum dot and the matrix. Therefore the Hamiltonian 
matrix h can loose hermiticity unless the proper recipe for the order of differential and 
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multiplication operators is chosen. This recipe is given by 

1 



f(r)didj 



fir)* - - (fit/(r) + /(r)fit,) . 



(9) 



All the elements in the Hamiltonian matrix are a linear combination of the elements of 
the form E l = f(r)eij, E 2 = f(r)e i:j ki, E 3 = f(r)x a y p , E A = f(r)x a y (3 k i and E 5 = 
f(r)x a y t3 k i kj, where k{ (i G {1,2,3}) is the differential operator ki = e-ij are the 

components of the strain tensor and f(r) is of the form 

f(r) = f QD X QD (r) + f M (l - XQD (r)), (10) 

where f® D is the value of a material parameter in the quantum dot and f M its value in the 
matrix, Xqd{t) is the quantum dot characteristic function equal to 1 inside the dot and 
outside the dot. Their Fourier transforms (jSJ) are within the recipe ([9]) thus given by 
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where Af = f M — f® D , eij{q) are Fourier transforms of the strain components given by 



Xqd (q) i s Fourier transform of the quantum dot characteristic functions 

x x qd (q) = 7^3i J QD ^ r & V ex P ' r ) , 

where the integration goes only over the volume of the quantum dot and 
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where a and (3 are non-negative integers from the set 

(a, (3) G {(0, 0); (1, 0); (2, 0); (0, 1), (1, 1), (0, 2)}. (19) 

The center of the pyramid base is taken as the origin of the coordinate system (Fig. [1]). The 
analytical formulae from which e^q) can be derived in a crystal with cubic symmetry are 



given in Ref. |33j. After integration, the characteristic functions can all be expressed as a 



linear combination of integrals of the type 



I m (q) = / x m e iqx dx, (20) 



o 



where m G {0, 1, 2, 3} and can be therefore evaluated analytically. All the integrals J x a y {i(q — 
k) are evaluated analytically, as well. 

In order to find the energy levels and the wave functions in the quantum dot, the eigen- 
value problem (CO) should be solved. The direct application of this approach would lead to 
an eigenvalue problem of a matrix of size N x N. However, it is possible to exploit the 
symmetry of the model to block diagonalize the corresponding matrix, as will be done in 
SecUD 



II. THE SYMMETRY OF THE MODEL 

As already discussed in the introduction, the symmetry group of the model is the double 
group C4. The generator of the group is the total angular momentum F z and therefore 
the representations of the elements of the group are given by the operators D(R(<p)) = 
exp(—ipF z ), where <p G {kir/2} {k G {0,1,..., 7}) and R(ip) is a rotation by an angle 
(p. In order to find how operators D(R(<p)) act on the states (Eq. it is enough to find 
how D(R(n/2)) acts on the states since D(R(kn/2)) = D(R(n/2)) k . The total angular 
momentum F z is a sum of the total angular momentum of the Bloch function J z and the 
orbital angular momentum of the envelope function L z , i.e. F z = J z + L^-. Therefore action 
of the operator D(R(tt/2)) on state |^) is composed of a rotation of the envelope functions 
in real space generated by its orbital angular momentum L z and a rotation of the Bloch 
function generated by its total angular momentum J z 

8 

D(R(tt/2))\V) = £ [exp(-z^4)^(^)] [exp(-i<pJ s )\i)] . (21) 
i=i 

Since the basis of Bloch states \i) is the eigenbasis of J z it follows that 

exp(-ij z ii/2)\i) = exp(-iJ z (i)n/2)\i), (22) 

where J z (i) is the eigenvalue of the z— component of the total angular momentum of Bloch 
function |i> (J t (l) = -1/2, J z {2) = 1/2, J z {3) = -3/2, J Z {A) = -1/2, J g (5) = 1/2, 
J z {6) = 3/2, J z (7) = -1/2, J z (8) = 1/2). Thus the operator D{R(tt/2)) acts on the state 
|^) as 

D( J R( 7 r/2))|*) = X;^(3/ > -^^)exp(-iJ,(i)7r/2)|i). (23) 
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By assuming the envelope functions as a linear combination of a finite number of plane 
waves, we have already reduced the otherwise infinite Hilbert space of the model to the 
Hilbert space Tt of dimension N formed by linear combination of plane waves multiplied by 
the Bloch functions. The basis of the space Ti is given by 

\k,i) = exp(ik ■ r)\i), (24) 

where k = 2tx (j^, j^, ffj, m x G {-n X} . . . , n x }, m y G {-%, • • • , n y }, m z G {-n z , ...,n z } 
and i G {1, 2, . . . ,8}. 

We shall first represent the operator D(R(tt/2)) in the plane wave basis of the space Ti.. 
We thus need to know how D(R(7r/2)) acts on the basis vectors. Using (1231 and ff24l) we 
find that 

D(R(n/2))\(k x ,k y ,k z ),i) = exp(-iJ z (i)n/2)\(-ky,k x ,k z ),i). (25) 

One can note that for (k x , k y ) = (0, 0) the acting of D(R(n/2)) on the basis vector is just a 
phase shift and the orbit of acting of the group elements is just an one-dimensional space (we 
shall denote it as W(o,o,fc*),i)> while for (k x , k y ) ^ (0, 0) the orbit is a four- dimensional space 
(that we shall denote as Ti.(k x ,k y ,k z ),ij where k x > and k y > to avoid multiple counting of 
the same space) with the basis 

|61) = \(k x ,k y ,k z ),i), 
\b2) = |( ky,k x ,k z ),i), 
1 63) = \(-k x ,-k y ,k z ),i), 

|64) = \(k y ,-k x ,k z ),i). (26) 

In the space 7^(o,o,fcAi the representation D reduces to an one-dimensional representation 
defined by 

D (Oi o, fc .),i0R(7r/2)) = exp(-iJ z {t)7t/2), (27) 

while in the space 7~C(k x ,k y ,k z ),i it reduces to a four-dimensional representation which is given 
in the basis from Eq. [26] by 



■ D (* B ,* w ,fc,),t(-R('r/2)) = exp(-iJ 2 (z)7r/2) 



Since the spaces / H(k x ,k y ,k z ),i and W(o,o,fc a ),i are invariant for the representation D, it is given 
by an orthogonal sum 

^=0 A^«,i + ©Aw,i- ( 29 ) 

From Eqs. [57|and[2H]one finds that the characters of the representation of the group elements 
are given by 

| 4 k = 

X (A^,M,^(W2))) = -4 k = 4 (30) 

[ fe G {1,2,3,5,6,7} 



1 
10 
10 
10 



(28) 
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and 

X (Ao,o,*.),i(fl(W2))) = exp(-iJ z (i)W2)- (31) 

The characters of the irreducible representations of the double group C4 are given by 
x(A mf (R(kir/2))) = exp(ikm f ir/2), where m f G {-3/2,-1,-1/2,0,1/2,1,3/2,2} and 
fc G {0, 1, ... , 7}. One finds from (EDI) that 

D(k x ,k y ,k z ),i = A x /2 + ^-1/2 + ^3/2 + ^-3/2 (32) 

and from (I3T]) obviously 

D(o,o,k x ),i = A-.j z (i). (33) 
Using ([32]) and ([33D, it follows from ([29]) that 

/) = + iViA_ 1/2 + N 2 A 3/2 + N 2 A_ 3/2 , (34) 

where 

Ni = 8n x (n y + l)(2n z + 1) + 3(2n z + 1) (35) 

and 

N 2 = 8n x (n y + l)(2n 2 + 1) + 2n z + 1. (36) 

Projection operators^ were then used to find the symmetry adapted basis. The projection 
operators are given by 

P Amf ((k x ,k y ,k z ),i) ^ExK(%/2)p ( WJ,i(%/2)), (37) 

8 fc=0 

while 

P Amf ((0,0,k z ),i) = 1 (38) 

and they project arbitrary states in space Ti to the elements of the symmetry adapted basis. 
The explicit forms of the projection operators are derived from fl37|) and (128]) and in the 
basis fT26|) are equal to 
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(39) 

where (k x ,k y ,k z ) was omitted in all brackets for brevity and where the matrices Mi, M 2 , 
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M 3 and M 4 are given by 
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The elements of the symmetry adapted basis are finally given as: 

|(0,0, W> ie {1,4,7} 

|(|61> + |62) + |63> + |M» ie {1,4,7} 

i(|61) -i\b2) - |63)+i|64)) ^ G {2,5,8} 

§(|W) - \b2) + \b3) - |M» 2 = 6 

i(|61)+i|62)-|63)-i|64)) i = 3, 



|A 1/2 ,(0,0,y,i) 

1-^1/2) (&xj ^z), 

1-^-1/2) {k x , k y , k z ), i) 

1-^1/2) (&x, fcj/j ^z), 

|-4i/2, {k x , k y , k z ),i) 



(40) 



(41) 



1^-1/2,(0,0 

1^.-1/2, {k x , ky 
1-4-1/2, {k x , k y 
1-4-1/2, {k x , k y 
1-4-1/2, {k x , k y 

|-4- 3/2 ,(0,0 

1-4-3/2, {k x , ky 
1-4-3/2, {k x , ky 
1-4-3/2, {k x , ky 
1-4-3/2, {k x , ky 



k z ),i) 
k z ),i) 
k z ),i) 
k z ),i) 
k z ),i) 

k z ),i) 
k z ),i) 
kz),i) 
k z ),i) 
k z ),i) 



1(0,0,*;,),;) i G {2, 5, 8} 

i (|61> + |62> + |63> + |M» 2G{2,5,8} 

i(|61)-i|62)-|63)+i|64)) 2 = 6 

|(|61>-|62> + |63>-|M» 2 = 3 

I(|61)+i|62>-|63>-<|64» <G {1,4,7}, 

|(0, 0, fc z ),i) 2 = 6 

i (|61> + |62> + |63> + |M» ? = 6 

i(|61)-i|62)-|63)+i|64)) 2 = 3 

|(|61>-|62> + |63>-|M» 2 G {1,4,7} 

| (|61) + i\b2) - |63> - i|64)) 2 g {2, 5, 8}, 



(42) 



(43) 



|-4 3 /2, (0,0, k z ), 2) 
1-43/2, (&x, fcj/j k z ),i) 
|-4 3 /2, {k x ,k y ,k z ),i} = \ (|61) 

1-43/2, (&x, fcj/j ^z), 

1-43/2, (&x, ^2/, k z ),i) 



|(0,0,^), i) 
i(|61> + |62> + |63> + |M» 



2' = 3 
2' = 3 

, i|62) - |63)+i|64)) 2 G {1,4,7} 

|(|61> - |62> + |63> - |M» 2G{2,5,8} 
i(|61)+i|62)-|63)-i|64)) 2 = 6. 



(44) 



The Hamiltonian matrix elements between basis elements having different symmetry are 
equal to zero implying that in this basis the Hamiltonian matrix is block diagonal with four 
blocks of sizes A^i x JV b ^ x ^, JV 2 x N 2 and N 2 x N 2 , respectively. The time necessary 
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to diagonalize the matrix of size N x N scales approximately as N 3 . Therefore the block 
diagonalization obtained reduces the computational time approximately by a factor of 16. 
Since all the basis vectors of the symmetry adapted basis are linear combinations of one 
or four vectors of the plane wave basis, it follows that Hamiltonian matrix elements in the 
symmetry adapted basis can be expressed as linear combinations of one, four or sixteen 
Hamiltonian matrix elements in the plane wave basis. Therefore all the elements of the 
four blocks are given by analytical formulae that can be easily derived from the analytical 
formulae for the matrix elements in the plane wave basis given in Sec. [H 

It can be proved by considering the two dimensional irreducible representations of the 
more general double group that states with the same absolute value of m/ are de- 

generate in pairs (the Kramer's degeneracy). This degeneracy is lifted in the presence of 
external axial magnetic field B when the symmetry reduces from C± v to C4 and the time 
reversal symmetry relation E mf (B) = E^ mf (—B) holds then. 

It has been pointed out 3 ^ 1 ^ that piezoelectric effects in single dots of realistic sizes are 
small, changing the eigen-energies of the system by less than 1 meV, and can be neglected. 
It has also been shown that in a vertically stacked double quantum dot the influence of 
a piezoelectric field is more important since the piezoelectric potential generated by the 
two dots adds up in the regions above and below the dots, while it is almost cancelled 
out in the region between the dots. Consequently, it is expected that in a periodic array of 
vertically stacked quantum dots considered in Sec. IHII the piezoelectric potential in the region 
between the dots would be almost cancelled out and that the influence of piezoelectric effect 
on eigen-energies would be small. Therefore the small piezoelectric potential that breaks the 
symmetry of the system from C4 to C2 can be treated as a perturbation. It belongs to the 
A 2 representation of C4 group and therefore only the matrix elements between the states 
with Arrif = 2 are non zero. Consequently, the piezoelectric potential doesn't change the 
energies in the first order of the perturbation theory and second order perturbation theory 
is needed to take the piezoelectric effect into account. 

We introduce the following notation for the electron states ne mf , where n is a positive 
integer labeling the states with given irtf in increasing order of their energy. Since the 
quantum number rrif originates from the irreducible representations of the double group C4 
whose elements are generated by the total angular momentum, its physical interpretation 
is that it represents the total quasi-angular momentum. The hole states will be labeled by 
nh mf , with the same meaning of the symbol as for the electron case, except that n labels 
the states in decreasing order of their energy, as is natural for holes. 

The Hamiltonian of the interaction with the electromagnetic field is obtained by replacing 
k with k + |A in the kinetic part of the Hamiltonian. In the dipole approximation A can 
be considered constant in space and furthermore all the terms quadratic in A are neglected. 
After calculating the matrix elements between states with a well defined symmetry, we 
obtain the following selection rules. If the light is z— polarized then Arrif = 0, while if the 
light is a^ 1 circularly polarized then Amj = ±1 (where by definition 3/2 + 1 = —3/2 and 
-3/2-1 = 3/2). 

III. RESULTS 

The method presented was applied to the calculation of the electronic structure of periodic 
array of vertically stacked pyramidal self-assembled quantum dots (Fig. [2]). The dimensions 
of the dots in an array were taken to be equal to those estimated for the structure reported 
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in Ref. [l9| - the base width b = 18 nm, the height h = 4nm, the wetting layer width 
dwL — 1.7ML, while the period of the structure in z— direction L z was varied in the interval 
from L z = h + dy/i where the dots lie on top of one another to L z = 16 nm. The dimensions 
of the_embedding box L x = L y = 2b were taken. The material parameters were taken from 



Ref. |36|. According to Bloch's theorem, the k— th component of the state spinor is given by 

tffc(r) =exp(iK z z)ip k (r), (45) 

where k G {1, 2, . . . , 8} and ipk(f) is periodic in the z— direction with the period L z . There- 
fore, the matrix elements in the four blocks of the Hamiltonian matrix for calculating E(K Z ) 
are given by linear combinations of the elements obtained by the same formulae from Sec. [I] 
except that k z and q z should be replaced by k z + K z and q z + K z , respectively. Since the 
relation E(K Z ) = E(—K z ) holds, only the states with K z > will be considered. The InAs 
unstrained conduction band edge is taken as the energy reference level. 




WL 



FIG. 2: Geometry of a periodic array of pyramidal quantum dots. The width of the pyramid base 
is b, the height h, the wetting layer width dwL, the period of the structure is L z . 

The small piezoelectric effect was assumed to be negligible in the calculation. In order to 
check this assumption it was included in the framework of second order perturbation theory 
(Sec. |TT])- Its influence on the state energies was of the order of lmeV and less, confirming 
our assumption. 

Two main factors influence the electronic structure of the periodic array of quantum 
dots - strain distribution and quantum mechanical coupling. 

The influence of quantum mechanical coupling is intuitively clear - as the distance between 
the dots increases the coupling is weaker implying smaller miniband widths. Due to their 
large effective mass, heavy-holes are the least influenced by coupling and the minibands of 
dominantly heavy-hole like states are narrow, while the minibands of electron and light-hole 
states are much wider. 

On the other hand, the strain distribution is complex and in principle all six components of 
the strain tensor influence the electronic structure. Still, the most important are hydrostatic 
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strain = en + e 2 2 + e 33 that determines the position of the electron and hole levels and 
biaxial strain = — ~ (en + 622) whose main influence is on splitting of the light and 
heavy-hole states. The bigger the value of hydrostatic strain, the lower the conduction band 
states are in energy and the higher the valence band states are in energy. When the biaxial 
strain is negative, the light-holes tend to have higher energy than heavy-holes, while when 
it is positive the situation is opposite. Having the importance of those two components of 
strain in mind, we have investigated first how they change when the distance between the 
dots in an array is varied. We have found that as the distance between the dots increases, 
the hydrostatic strain in the dots decreases. On the other hand, for small values of the 
period, the biaxial strain is negative, while for larger values it changes sign and increases 
further. 



A. Energy levels in the conduction band 

The dependence of the miniband minima and maxima on the period of the structure 
is given in Fig. [3j This behavior is expected. When the dots are close, they are strongly 
coupled and the minibands are wide, while as L z increases the coupling is weaker and the 
energy spectrum becomes discrete. For large values of the period when the miniband width 
practically vanishes, we still see a rising trend in energy. This rise is caused by a still 
decreasing value of hydrostatic strain. This leads us to the conclusion that the range of 
strain effects is larger than the range of quantum mechanical coupling. 
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FIG. 3: The conduction miniband minima and maxima dependence on the period of the structure 
L z , The le±i/2 miniband is represented by squares, the le±3/2 and 2e±i/2 miniband by circles, 
the 2e±3/2 miniband by triangles. 

The ground miniband has \mf\ = 1/2 symmetry, while the first and second excited 
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miniband having different symmetries \m,f\ = 1/2 and \m,f\ = 3/2 are nearly degenerate. 
Their difference in energy is less than 1 meV, too small to be seen on the graph. A comment 
should be given about the near degeneracy of those two states. It has been practice in the 
literature to say that these two states are exactly degenerate in the absence of a piezoelectric 
effect and that the piezoelectric effect breaks the degeneracy of these states. This is indeed 
true if the carrier energy spectrum in the quantum dot is modeled by simple one band 
Schrddinger equation. The symmetry group is then C^ v , consisting of transformations in 
real space generated by the orbital angular momentum L z and the first and second excited 
state transform according to the same two dimensional irreducible representation of C^ v 
implying their degeneracy. However, when the 8-band k ■ p model, which is inherently spin- 
dependent, is used, the Hamiltonian no longer commutes with the rotations generated by 
orbital angular momentum, but the total angular momentum. The symmetry group is the 
double C4 group that has only the one dimensional irreducible representations and there is 
no a priori reason for the states with different absolute values of mj to be degenerate. 

In order to explain the near degeneracy of the states with different symmetry, we have 
further investigated the 8-band k ■ p Hamiltonian and checked that it would still commute 
with the transformations in real space generated by orbital angular momentum if the valence 
band spin-orbit splitting would be set to zero. Since the influence of valence band spin-orbit 
splitting on the levels in the conduction band is not substantial, the degeneracy of the 
two states is small. We thus conclude that the origin of splitting of the first and the second 
excited state is not just the piezoelectric effect but also the valence band spin-orbit splitting. 

All the minibands shown in Fig. [3] exhibit minima at K z = and maxima at K z = 7t/L z 
for all values of the period L z . For small values of the period L z there is an energy overlap 
between different minibands, while the minibands are completely separated for larger values 
of L z . There is no crossing between states of different symmetry. 

B. Energy levels in the valence band 

The miniband minima and maxima dependence on the period of the structure are given 
in Fig. H] for the three highest minibands in the valence band. Due to the combined effects 
of strain, mixing of light and heavy-holes and quantum mechanical coupling between the 
dots, the hole minibands exhibit a more complex structure than the electron minibands. 

In order to explain such behavior we note first that the effective potential felt by carriers 
depends on K z . We define the effective potential as the value obtained by diagonalizing the 
Hamiltonian with k x = k y = and k z = K z . The z-dependence of electron, light, heavy and 
spin-orbit split hole K z = and K z = ir/L z effective potentials at the pyramid axis for a 
few different values of L z is shown in Fig. 

Since the effective mass of the light-holes is small, the light-hole effective potential is 
substantially different for K z = and K z = tt/L z , while in the case of heavy- holes that 
difference is much smaller. As seen from Fig. [5j as the period of the structure increases, the 
effective potential felt by K z = light-holes decreases, while quite oppositely the effective 
potential felt by heavy-holes increases. Both of these trends are an expected consequence of 
the increase in the value of the biaxial strain. Consequently, in the range of low values of 
L z the hole states with K z = are dominantly of the light-hole type, while the states with 
K z = tt/L z are dominantly heavy-hole like. The states with K z = tt/L z remain of heavy- 
hole type across the whole investigated interval of L z and their energy therefore increases 
with increasing L z . The energy of the light-hole K z = states decreases with increasing L z 
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FIG. 4: The miniband minima and maxima dependence on the period of the structure L z . The 
l/i±i/2 miniband is represented by squares, the l/i±3/2 miniband by circles and the 2h± 3 / 2 miniband 
by triangles. 
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FIG. 5: Effective potentials at three different values of the period L z at K z = (full lines) and 
K z = ir/L z (dashed lines) for electrons (E), light-holes (LH), heavy-holes (HH) and spin orbit split 
holes (SO). 

and at the same time their heavy-hole content increases. As a consequence of two different 
energy trends for K z = and K z = n/L z the miniband width decreases for all the states 
until a certain point where the energy of the K z = state becomes less than the energy of 
the K z = tt/L z state. This point, where the inversion of the sign of the miniband effective 
mass occurs is different for different states, for the l/i±i/2 state it occurs around L z = 6.5 nm, 
while for the lh± 3 / 2 and 2/i ±3 / 2 states it occurs around L z = 5.8 nm. The light-hole content 
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of the K z = states decreases with increasing L z and eventually they become dominantly 
heavy-hole like. This light to heavy-hole transition occurs at L z = 7.5 nm for the l/t±i/2 
state, at L z = 5.5 nm for the l/i±3/2 state and L z = 5.0 nm for the 2h± 3 / 2 state. Since the 
energy of the light-hole states decreases with increasing L z and the energy of the heavy-hole 
states increases, the position of the light to heavy- hole transition corresponds approximately 
to the position of the energy minima of K z = states, as can be verified from Fig. HI We 
further observe that the ground hole state for L z < 6.5 nm is l/t±i/2 having |m/| = 1/2 
symmetry, while for L z > 6.5 nm it is l/i±3/2 having |m/| = 3/2 symmetry. Therefore at 
the critical point L z = 6.5 nm, we observe an interesting effect of a simultaneous change of 
ground hole state symmetry, a change of the sign of the effective mass and a change of the 
ground state type from light to heavy-hole like. 

The spin orbit split band certainly influences the exact positions of the energy levels, 
however, being far in energy from the light and heavy-hole bands it doesn't influence the 
overall behavior described in the previous paragraph. This is verified by the fact that the 
spin orbit band content of the hole states is typically of the order of 5%. 

As far as the spatial localization of the wave functions is concerned, one would expect 
from the effective potential profiles given in Fig. [5] that dominantly light-hole states would 
be confined outside the dots and dominantly heavy-hole states inside the dots. However, the 
states are of light-hole type only when the dots are very close to each other and the effective 
potential well is then too narrow to confine the hole. Therefore, the light-hole like states are 
spread both inside and outside the dots. When the distance between the dots increases and 
light to heavy-hole transitions take place, the hole state becomes localized inside the dots. 

C. Influence of external axial magnetic field 

The magnetic field dependences of the miniband minima and maxima of the conduction 
and valence band states for the structure with the period L z = 6nm are shown in Figs. [6] 
and El respectively. As already mentioned, Kramer's degeneracy is broken in magnetic field. 
The relation E mf (B) = E_ mf (—B) holds, thus only the B > part of the dependence is 
shown on the graphs. 

The magnetic field splitting between le+3/2 and le_3/2 states, as well as between 2e_i/2 
and 2e + i/2 states is significant because the mesoscopic angular momentum^ 3 , of those states 
is different from zero. However, the splitting between le_i/2 and le + i/2 states and between 
2e + 3/2 and 2e_3/2 is much smaller, too small to be seen on the graph (of the order of 
few meV). There is no crossing between the states of different symmetry. We also note 
that the le + i/2 and le_!/ 2 minibands overlap with 2e_i/2 and le_3/2 minibands for B = 
but as the magnetic field is increased this overlap vanishes (for B > 12 T). The energy 
separation of the minibands has an important affect on the dynamical characteristics of the 
structure since it suppresses all the one particle energy conserving scattering mechanisms 
between those minibands (like ionized impurity scattering) and with further separation even 
suppresses the mechanisms with energy exchange (like acoustic phonon and longitudinal 
optical phonon scattering). 

The splitting between the hole states is also of the order of a few meV, however since 
the energy difference between different states is also small, this splitting is enough to cause 
crossings between states of different symmetry. We further find from Fig. [7] that the mini- 
bands l/i_3/2 and I/1+3/2 that are degenerate at B = become completely separated already 
at B > 3T and the same effect for 2/i_ 3 / 2 and 2h +3 / 2 occurs at B > 23 T. Although the 
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B[T] 

FIG. 6: Magnetic field dependence of miniband minima and maxima for le_ 3 /2, le_i/2, le + i/2? 
le +3 / 2 , 2e_ 3 /2, 2e_ 1/ / 2 i 2e +1 / 2 and 2e +3 / 2 states. The position where different minibands separate 
is marked. 

magnetic field splitting is the most pronounced for lh_i/ 2 and lh + i/ 2 states, the effect of 
the separation of those minibands occurs at magnetic fields larger than 30 T, which is a con- 
sequence of larger miniband width than in the previous cases. Apart from separation of the 
minibands, the magnetic field can also concatenate otherwise nonoverlapping minibands. It 
is seen in Fig. [7] that \h + i/ 2 and l/i_3/2 start to overlap at B ~ 11 T and that for B > 22 T 
the range of energies of 1/2.-3/2 becomes a subset of the range of energies of l/i+i/2 miniband. 

IV. CONCLUSION 

In conclusion, we have developed a symmetry-based method for calculation of electronic 
states in pyramidal InAs/GaAs quantum dots. The corresponding Hamiltonian matrix ob- 
tained by the plane wave method was block diagonalized into four matrices of approximately 
equal size, which enabled significantly faster calculation of energy levels within the plane 
wave method. The symmetry considerations not only enabled more efficient calculation of 
the electronic structure but also give more insight about the physics of the model by intro- 
ducing the quantum number of total quasi-angular momentum and giving the selection rules 
for interaction with electromagnetic field. The method developed was applied to calculate 
the electronic structure of a periodic array of vertically stacked pyramidal self-assembled 
quantum dots. It was found that as the distance between the dots is increased, at a certain 
critical point the ground hole state simultaneously changes symmetry from \rrif\ = 1/2 to 
\m,f \ = 3/2 and type from light to heavy-hole. The influence of magnetic field on the energy 
levels is in general less pronounced than the influence of quantum mechanical coupling and 
strain but nevertheless it can be used for fine tuning of the properties of the structure since 
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FIG. 7: Magnetic field dependence of miniband minima and maxima for l/i_ 3 / 2 (full circles), 
1^+3/2 (empty circles), lh_i/ 2 (full squares), lh + i/ 2 (empty squares), 2/i_ 3 / 2 (full triangles) and 
2e_|_3/ 2 (empty triangles) states. The positions where different minibands separate or concatenate 
are marked. 



its increase or decrease is able to separate energy overlapping minibands. 
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